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Abstract. We describe a new analysis of Y(nS ) — > n + n~ decays collected in pp collisions with the CDF II 
detector at the Fermilab Tevatron. This analysis measures the angular distributions of the final state muons in the 
Y rest frame, providing new information about Y production polarization. We find the angular distributions to 
be nearly isotropic up to Y p T of 40 GeV/c, consistent with previous measurements by CDF, but inconsistent 
with results obtained by the DO experiment. The results are compared with recent NLO calculations based on 
color-singlet matrix elements and non-relativistic QCD with color-octet matrix elements. 



1 Introduction 

A recent analysis [ 1 1 of T — > p + pT decays collected with 
the CDF II detector provides new measurements of the dis- 
tributions of decay angles, which depend on the polariza- 
tion of Y states produced in p~p collisions. Previous mea- 
surements 112131 carried out by the CDF and DO experi- 
ments provided useful, but incomplete information about 
these angular distributions and did not strongly favor the 
predictions of any of the models 15 141 used to calculate 
the production cross sections. Furthermore, the fact that 
these previous measurements are in apparent disagreement 
has led to the speculation that significant acceptance bi- 
ases could have been overlooked, motivating the need to 
perform additional tests of internal consistency in future 
measurements |[6) . 

The analysis described here measures the full angular 
distributions of the final state muons from T(IS), Y(2S), 
and T(3S) decays as functions of the Y transverse mo- 
mentum up to 40 GeV/c. This is the first analysis to re- 
port measurements of the Y(3S ) spin alignment. It is also 
the first analysis to measure spin alignment in two differ- 
ent coordinate frames and to compare rotationally invariant 
quantities in these frames to demonstrate internal consis- 
tency. 



2 Analysis Overview 

In the rest frame of the Y decay, the direction of the posi- 
tive muon is described using polar angles (8, tp) measured 
with respect to a given set of coordinate axes. The s-channel 
helicity frame, used in earlier analyses, defined the z-axis 
along the Y momentum vector, with the x-axis in the pro- 
duction plane and the j/-axis perpendicular to both x- and 
z-axes. Alternatively, the Collins-Soper frame [7| can be 
used for which the z-axis approximates, on average, the di- 
rection of the velocity of the colliding partons. The form 
of the angular distribution is constrained by angular mo- 
mentum conservation; for a vector meson decaying to two 



fermions it can be written 



dN_ 

do 



1 + Ag cos 2 8 + A v sin 2 8 cos 2tp + Ag^ sin 28 cos tp + 



Ah sin 2 6 sin 2<p + Ag sin 26 sin tp. 



(D 



e-mail: mjones@physics .purdue . edu 



A four-fold symmetry in the acceptance may be exploited 
to increase statistics in small bins of solid angle by combin- 
ing (8, tp) with (8, -tp) and (8, tp) with (n-6, n-tp), although 
this leads to a cancellation of the terms with A^ and A^ as 
coefficients. The coefficients quantify the shape of the an- 
gular distribution and provide direct information about the 
polarization of the ensemble of Y states since they are re- 
lated to the elements of the spin density matrix elements 
pi) AS) by the expressions Ag = (p n - poo)/(pn + Poo), 
=Pw/(Pu +Poo), and A eifJ = p 1 _i/(pn +p o)- 
When a sample of Y decays is selected using a dimuon 
trigger, the observed angular distribution will differ from 
the form in Eq. ([TJ because of the limited acceptance im- 
posed by the muon pj thresholds in the trigger and the ge- 
ometric coverage of the detector systems. The acceptance 
can change rapidly with both the transverse momentum 
and mass of the dimuon system but can be calculated ac- 
curately using a combination of Monte Carlo simulations, 
which model the detector geometry, and trigger efficien- 
cies measured using independent data samples. The an- 
gular distributions of dimuons with mass that include one 
of the Y(nS ) resonances will depend strongly on angular 
distributions present in non-resonant backgrounds, which 
can be highly non-isotropic. These may have angular dis- 
tributions that can be parameterized using Eq. ([TJ, but they 
could be more complex since they do not necessarily arise 
from the decay of a single vector state. 



2.1 Previous Analyses 

Previous Y spin alignment analyses I2I3I were performed 
using only the s-channel helicity frame and integrated an- 
gular distributions over tp, retaining only sensitivity to the 
coefficient Ag which is frequently denoted a in the litera- 
ture. This was carried out in several ranges of pr(Y) by 
fitting the dimuon mass distribution in discrete ranges of 
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cos 8 to determine the Y yields, correcting for detector ac- 
ceptance, and fitting the resulting distributions to a func- 
tion of the form 1 + a cos 2 8. In practice, this is not a trivial 
procedure because both the T acceptance and the shape of 
the background mass distribution change significantly with 
both pr and cos 8. 

The limitations of these analyses have been pointed out 
for several years now |6|. The measurement of only one of 
the three coefficients in Eq. ([TJ does not allow the calcu- 
lation of rotationally invariant quantities or the transfor- 
mation Ag from the s-channel helicity frame to different 
coordinate systems. This procedure also does not general- 
ize well to the analysis in many small bins of cos# and 
tp because the large number of fits to invariant mass dis- 
tributions with poorly constrained background shapes may 
suffer from large statistical fluctuations or systematic bi- 
ases. 



2.2 A New Approach 

The analysis procedure described here divides (cos 8, tp) 
into 20 x 36 bins, but it avoids the need to measure T yields 
in each bin separately. Instead, all dimuon events with in- 
variant mass near each of the T signals are selected and the 
observed numbers of events in each bin are modeled using 
separate angular distributions of the form in Eq. (Q]i for sig- 
nal and background, multiplied by the detector acceptance 
which is calculated in each individual bin. The parame- 
ters Ag, Ay, Ag^ corresponding to the angular distribution of 
the T signal can then be measured provided the amount of 
background and its angular distribution are known. 

Within a given range of dimuon pj, the amount of 
background under the T(nS) signals is determined from 
a fit to their invariant mass distribution but an independent 
sample is needed to constrain the shape of the angular dis- 
tribution from background sources. Such a sample is ob- 
tained by demanding that the extrapolated trajectory of at 
least one of the muons misses the average beam axis by a 
distance |o?ol > 150 pm. Although this "displaced" muon 
sample contains a few percent of the T signal due to the do 
measurement resolution, it mostly selects muons produced 
in semileptonic decays of heavy quarks, which forms the 
dominant source of background. Since the impact param- 
eter requirement does not bias the muon decay angle, the 
angular distribution of muons from background sources is 
expected to be the same in the complementary "prompt" 
muon sample. 

A simultaneous fit is then performed to the angular dis- 
tributions of dimuon events in prompt and displaced sam- 
ples selected from ranges of invariant mass around each of 
the T signals. In each bin of Pt(Y), the fit is performed by 
maximizing the likelihood function constructed from the 
probabilities of obtaining the observed numbers of events 
in bins of (cos 8, tp) given expected yields calculated using 

dN„ 

— N r fp$lr(cos 0,-, <pj) ■ w r (cos 6 h <pj; A r ) + 

CtLJ[j 

N d s p Jl b (cos 8 U <pj) ■ w h {cos 8i, tp } \ A b ), (2) 

dN d 

— N r (l - fp)Mr(cos 8,, tpj) ■ w T (cos 8(, tpf, Ay) + 

CtLJ[j 

N d $i b (cos 8i, (pj) ■ w b (cos 6 U <pf, A b ). (3) 



In these expressions, Nr and Nd are the numbers of T and 
displaced background events, f p is the fraction of the T 
signal retained in the prompt sample, and s p is the ratio 
of the background yields in prompt and displaced sam- 
ples. The parameters f p and s p are constrained using fits 
to the prompt and displaced mass distributions. The ac- 
ceptance J\y for T signal and Ji b for dimuon background 
events are calculated using Monte Carlo simulations and 
the measured trigger and muon selection efficiencies. The 
underlying angular distribution of muons from T decays, 
wj, is calculated using Eq. (JTJ with coefficients denoted 
collectively as Ar- The angular distribution of muons in 
the background component, w b , is similar, but has an addi- 
tional term, described below, that allows a better descrip- 
tion of the data. 



3 Analysis of CDF Data 
3.1 Upsilon Trigger 

The analysis of T — > p + p~ decays is performed at CDF us- 
ing a sample of events collected with a 3-level dimuon trig- 
ger. This trigger required the presence of two oppositely 
charged tracks at level 1 with pj > 1.5 GeV/c that extrap- 
olate to hits in one of the CDF muon detector systems [9 1. 
At least one of the muons had to be in the central region 
and the level 2 trigger required that it was also detected in 
a second muon detector system located behind additional 
steel absorber. After full event reconstruction, the level 3 
trigger required that this muon had pj > 4 GeV/c, the 
other muon had pj > 3 GeV/c and that the invariant mass 
of the pair was between 8 and 12 GeV/c 2 . The geometric 
acceptance of these triggers restricts the rapidity of the T 
sample to the central region, \y(T)\ < 0.6. 

For most of Run II, the level 2 trigger was prescaled 
dynamically to maintain an approximately constant accept 
rate and dead time, but more recently, the level 1 trigger 
was disabled when instantaneous luminosities were greater 
than 280 x 10 3() cm _2 s _1 . The prescaled triggers integrated 
approximately 70% of the delivered luminosity, averaged 
over the first 6.7 fb of data collected in Run II. Figure Q] 
shows the mass distribution of dimuons collected using 
these triggers that are used in the angular analysis. 



3.2 Angular Distributions in Background Events 

The analysis of the angular distributions of muons from T 
decays relies on accurately subtracting the angular distri- 
butions that are present in the background which are esti- 
mated using the displaced track sample. The validity of this 
procedure is checked by comparing the observed distribu- 
tion of decay angles in the high and low mass sideband 
regions. The example shown in Fig. [2] demonstrates that 
this is the case, with consistency between the two samples 
tested by computing the Kolmogorov-Smirnov statistic. 



3.3 Simultaneous Fit in Signal Regions 

The similarities observed in the angular distributions of 
prompt and displaced muon samples at masses both above 
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Fig. 1. Distribution of pfyT invariant mass. Pairs where both orig- 
inate from close to the beam axis are indicated by the solid his- 
togram, while those for which at least one misses the beamline 
by more than 150 yum are shown in the dashed histogram. The 
shaded regions indicate the range of masses used to select the 
Y(IS), T(25) and Y(3S) states, while the other regions are used 
to study background properties. 
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Fig. 2. Comparisons of cos (a,b) and <p (c,d) distributions, mea- 
sured in the Collins-Soper frame, for prompt and displaced sam- 
ples (histograms and error bars, respectively). Reasonable agree- 
ment, quantified by computing the Kolmogorov-Smirnov statis- 
tic, is observed in mass regions both below (a,c) the T(15) and 
above (b,d) the T(35) resonances. 
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Fig. 3. Projections of angular distributions measured in the 
Collins-Soper frame (a,b) and the .s-channel helicity frame (c,d) 
for muon pairs with 4 < p T < 6 GeV/c and mass in the vicin- 
ity of the T(1S) signal. The projected distributions in the data 
are indicated with error bars for the prompt (black) and displaced 
(gray) samples. Solid histograms show the projections of the fits. 



The quality of the resulting fit is assessed by comparing 
projections of the angular distributions observed in the data 
with the corresponding projections of the fit. Fig. [3] shows 
an example of projections for muon pairs selected from 
the mass range containing the Y(IS) signal. The quality 
of fits applied to the other kinematic regions is similar and 
shows no systematic trends that depend on either pr(p: + pr) 
or their invariant mass. 



3.4 Results 

The new measurements of Ag for the Y(IS) state can be 
compared with previous results from the CDF and DO ex- 
periments and with recent NLO predictions II 1 01 1 11121 . These 
are shown in Fig.|4]from which it is apparent that although 
the new results are consistent with the Run I CDF mea- 
surement, they are inconsistent with DO analysis, with the 
significance estimated to be approximately 4.5<x. The theo- 
retical predictions are currently somewhat imprecise due to 
the poorly measured production cross sections for \bj{nP) 
states that decay into Y(IS) [13], but recent results from 
the ATLAS experiment lTBl may improve on this situation. 



and below the Y(nS ) resonances support the use of the dis- 
placed muon sample to model the background properties 
under the Y signals. Thus, we apply the simultaneous fit 
to prompt and displaced muon samples with mass selected 
from each of the three regions containing the T signals in- 
dicated in Fig. ([1]). To accommodate the description of the 
angular distribution of the background, which is observed 
to be very non-isotropic, an additional term proportional to 
cos 4 is added to Eq. (Q} to obtain the functiona used for 
u>b- In addition, a component of the sample that is strongly 
peaked at large values of cos 9 in the s-channel helicity 
frame is removed by requiring that \pr{p- + ) - Pt(m ~)l < 
(pr(ji + fJ.~) - 0.5 GeV/c). This restriction is included in the 
calculation of the T acceptance but has a negligible effect 
for p T (Y) > 6 GeV/c. 



4 Rotational Invariants 

We also demonstrate the internal consistency of the results 
by calculating the rotational invariant A = (Ag + 3A ip )/(l - 
A^ using the values of Ag and A v measured in the s-channel 
helicity frame and the Collins-Soper frames. Agreement 
between the values calculated in each coordinate frame is 
an important consistency test because poor determination 
of the experimental acceptance or inaccuracies in the sub- 
traction of the highly non-isotropic backgrounds would be 
expected to introduce coordinate frame dependent biases 
in the measured angular distributions. 

The value of A quantifies the shape of the angular dis- 
tribution independent of its orientation with respect to a 
coordinate frame. Decays of T states with pure transverse 
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Fig. 4. Measurements of in the ^-channel helicity frame for 
the T(15) state compared with previous results from CDF (with 
\y(T)\ < 0.4) and DO (with \y\ < 1.8) and with next-to-leading 
order calculations based on NRQCD with color-octet matrix ele- 
ments, and a next-to-leading order color-singlet model. 
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Fig. 5. Rotational invariant A, measured in the Collins-Soper 
frame (dark lines) and .s-channel helicity frame (gray lines) as 
functions of pr(X). 



polarization yield angular distributions with A = +1 while 
a purely longitudinal polarization gives A = - 1 . A value 
of A = correspond to an isotropic angular distribution, 
which cannot result from the decay of a pure spin-1 state 
but instead would indicate that multiple production mecha- 
nisms are present leading to an effectively unpolarized en- 
semble of decays. 

Figure[5] shows the values of A measured for the Y(1S ), 
Y(2S) and Y(3S) states as functions of pr(Y) in both the 
Collins-Soper and the s-channel helicity frames. Since A 
is measured in each coordinate frame using the same data 
samples, the statistical uncertainty of each measurement 
is highly correlated. The sizes of variations in values of A 
measured in the two frames that would be expected from 
purely statistical fluctuations were estimated using Monte 
Carlo simulations and found to be generally consistent with 
the observed differences. This is the first analysis to per- 
form such a test and based on these findings, there appears 
to be no evidence for significant biases in the calculated 
acceptance. 



The values A » reached at large pj suggest that all 
three of the Y(nS) states are produced in an unpolarized 
mixture. This is the first measurement of angular distribu- 
tions in Y(3S) decays and is significant because it had been 
thought that a greater fraction of 3S states should be pro- 
duced directly, rather than via feed-down from Xb states, 
in which case the calculated spin alignment predictions 
should be more precise. 



5 Conclusions 

The measurements described here provide the most de- 
tailed characterization of the angular distributions of Y — > 
p + pT decays produced at a hadron collider to date. We find 
little evidence for strong polarization of any of the three 
Y(nS ) states in the central region of rapidity \y\ < 0.6 and 
with pj up to 40 GeV/c. This is consistent with the results 
previously obtained in Run I by CDF and inconsistent with 
measurements carried out by DO in Run II. Although the 
DO measurements were carried out over the wider range 
of rapidity \y(Y)\ < 1.8, we find no evidence that the an- 
gular distributions change rapidly in the central region of 
rapidity accessible to the CDF detector. We look forward to 
the possibility of refined predictions from theory and new 
results from the LHC experiments which may be able to 
further clarify the experimental situation. 
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